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FORCE  AT  A POINT  IN  THE  INTERIOR  OF  A SEMI-INFINITE  SOLID 


Introduction 

In  a paper ^ under  the  same  title,  the  solution  of  the  linear  equa- 
tions of  equilibrium  of  an  elastic  bod7  was  given  for  the  case  of  a force 
acting  at  a point  within  an  isotropic  body  bounded  by  a plane.  The  result 
was  obtained  by  starting  with  Kelvin’s  solution  for  a force  in  an  infinite 
body  and  guessing  the  nuclei  of  strain  to  add  outside  of  the  semi-infinite 
body  so  as  to  annul  the  tractions  on  the  plane  boundary.  In  the  present 
paper  it  is  shown  how  these  results  may  be  obtained,  directly,  by  means  of 
an  application  of  potential  theory. 

fopkovitch  Functions 


In  an  isotropic  elastic  body  in  equilibrium,  the  displacement  u is 
governed  by  the  equation ^ 

+ F =0  , (D 

M I A ]/  M ^ ' 

where  jjl  is  the  shear  modulus,  V is  Poisson’s  ratio  and  F is  the  body 
force  per  unit  of  volume. 

For  an  isotropic  body,  the  stress,  W , is  related  to  the  displace- 

ment  by 

CT  = /*V-uJ  + >u(Vu  t U V)  (2) 


^ R.  D.  Mindlin,  ifcyaics,  Vol.  7 (1936),  pp.  195-202. 

(2)  yor  the  vector  notation  used  in  this  paper,  see  C.  E.  Wsatherburn , 
Advanced  Vector  Analysis,  G.  Bell  and  Sons,  Ltd.,  London,  1928. 


where  A*  2//  J//(l -2  v). 

By  Helmholtz’s  theorem^  \ U may  be  resolved  into  lamellar  and  sole- 
noidal  components : 

u * , V-H  m0  } (3) 

so  that  (l)  may  be  written 

u vx(ct  Vif  + P*  w)  + F - 0 
where  ot  - 2 (\-vjj(\-2v)  . 


i.e. 


The  quantity  in  parenthesis  in  (4)  is  a vector,  say,  B 


r-(f  + V*  H = B , 

^ M IA  ^ / 


ytf  B * - F . 

^ VA  kA 

Operating  on  (5)  with  U-  , ve  find 

Ot  Vlv  «=  v*  B 

' ^ v» 

the  complete  solution  of  which  is 

2 oaf  « r-B 

M kA  1 


(4) 

iBx 

(5) 

(6) 

(7) 

(8) 


(9) 


where  ^ is  a scalar  function,  which  satisfies 

s * r f , 

and  f- lX  tjy  +kz  is  the  position  vector. 

•*  *»  v.  >A 

Substituting  (8)  in  (3)  and  eliminating  7*  H by  means  of  (5)»  there 

results 

Sf  “9"  + <*) 

^ 7aB  - -F 


(10) 

(11) 


(3) 


c/eatherbum,  p.  44. 
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uV'/S  = r-F  . (12> 

Thus  the  displacement  is  expressed  in  terms  of  the  Papkovitch  functions,  B 
and  ^ , whose  Laplacians  are  known  if  the  body  force  F is  known. 

The  proof  of  completeness  of  the  ftipkovitch  functions,  given  above, 
is  an  extension,  to  include  the  body  force,  of  one  given  in  a previous 
paper  ^ , where  0 and  @ were  called  Papkovitch'^  functions  after  the 
originator  of  the  solution  (10)  of  the  elasticity  equations.  Recent  writers 
associate  these  functions  with  the  name  of  Boussinesq^,  who  introduced 
and  0 , but  employed  functions  of  a different  type  where  B and  £3  could 

n j 

have  been  used. 

Green’s  Formula 

The  value  of  a function  V , at  any  point  in  a region,  may  be  ex- 
pressed in  terms  of  its  values  at  the  boundary,  its  Laplacian  and  Green's 

(7 ) 

function,  G , for  the  region,  by  means  of  Green’s  formula  1 

- +irV  = Jvn  VGdS  + J*G  V2Vdv  ^ 

For  the  region  z ^ 0 , Green’s  function  is 

G = r%~'  - r’'  (U) 

where 

^ R.  D.  Mindlin,  Bull.  Am.  Math.  Soc.,  Vol.  42  (1936),  pp.  373-376. 

P.  F.  Papkovitch,  Comptes  Rendus,  Acad,  des  Sciences,  Paris,  Vol.  195 
(1922),  pp.  513-515  and  754-756. 

^ J.  Boussinesq,  "Application  des  Potentials  a 1 'Etude  de  l'Equilibre  et 
du  Mouvement  des  Solides  Elastiques,"  Gauthier-Villars , Paris,  1885,  pp.  63 
and  72. 

(n\ 

' ' Weatherbnrn,  p.  34* 
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I i^r=. ■■■  II  ■ — — ■ i—— — I— — — — — — — — — — t— — ■ — — ■ I . I mm —— ww—  » 1 — 


cx  = (t-ff  + b-v)1  + (z-  s)x 


(15) 


r»a  - L*-s)x  +ki-jf  + fc  + *f  d*) 

in  which  (x,y,Z  ) are  the  coordinates  of  a point  ) in  the  region 

and  (f,  rr)>  5 )♦  5 ) are  coordinates  of  a source  point  ) 

and  its  image  <%'  y - $ ),  respectively. 


Force  at  a Point 

Kelvin's  definition  of  a force  at  a point  takes  the  following  foray 
in  the  present  case.  Consider  a distribution  of  body  forces  F in  a closed 
region  T within  Z > 0 , with  F * 0 outside  T but  within  z > 0.  uimirish 
T indefinitely!  always  enclosing  the  point  C (o,0}c  ),  but  let 

sb+ru  /*  F dv  = P (17) 

T-*o  - 

where  P is  a constant  force  at  C • 

For  later  use,  we  note  that  the  limit,  as  T approaches  zero,  of 


Qlfi,  "1,  ?)  » clo.»,c)  , 

(18) 

=C'(»,o,-c)  , 

(19) 

r1  ■=  x*  + v'  + (z-c)1  - R,J  ^ 

(20) 

15*  « x*  +y*+  (z  + t)1  * Rj  , 

(21) 

(22) 

9 / 1 1 \ 

Sz  (,  R + f?J  ’ 

(23) 

3a<x  __  3Z  1 1 l \ 

W ■ f?J  ’ 

(24) 

-4' 


I 


3 6 

3$ 


a / 1 i \ 


(25) 


Force  Normal  to  the  Plane  Boundary 

In  this  case  we  take  F = = 0 and  0X  = B^  - O . The  remaining 

Papkovitch  (in  this  case  Boussinesq)  functions*  B,  and  (9  * must  satisfy 
the  condition  of  vanishing  traction  on  z - O*  Thus*  we  have*  from  (2) 
and  (10)*  on  z = 0 


<r«  = 


W^l 


(26) 


(27) 


(28) 


The  function  in  brackets  in  (26)  is  one  whose  Laplacian 

_|if_ 1=- (29) 
|_l  ^32  3Z1  J M M 


is  known  throughout  Z ^ 0 and  whose  boundary  value  is  zero.  Hence*  by  (13 )* 


_ 3*/9  _ _l 


3z  7)Z 


TT  ~ 


47T/W 


. oo) 


Now*  integrating  the  first  term  in  the  volume  integral  by  parts* 

J]fG  If-  = JjGF.dM'? 
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_b 


The  surface  Integral  in  (31)  vanishes  because  =G  «0  on  the  boundary  of 
the  body.  Then*  by  (17)  and  (23 ), 


-fz~o  = p. 4z-(-^ 4 -^)  . 


(32) 


Similarly,  integrating  by  parts  twice  and  using  (17),  (18)  and  (24), 

f&  dV  = C **  "lzr  " T5T  ) * (33  ) 

Hence 

- If"  = Tn>r  [2M(-^  4 ^)  + c Srfc  ~ 

where  one  integration  with  respect  to Z has  been  performed.  (The  arbitrary 
function  of  x and  y , thereby  introduced,  must  vanish  since  2(\-v\&z~  Dpjdz 
must  vanish  as  Z— *<*»). 

Returning  to  the  boundary  conditions,  we  note  that  (27)  and  (28) 
can  be  integrated  with  respect  to  x fluid  y,  respectively,  so  that,  on  z ~ 0} 

(l-lv) Bz~  If-  =0  . (15) 

The  Laplacian 

[('- Jv)  B,  - ^ F,  - jr  ^ (16) 

is  known  throughout  z ^ O and  hence,  by  (13), 

(i-»te-3f-  = fv  ■ (,7) 

i I 
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By  the  same  process  as  before , using  (17),  (18),  (22)  and  (23),  we  find 

(l-2^)Bj  - If-  = R[)  + c.  * T^U  . (38) 

Finally,  from  (34)  and  (38), 


R _ S 
3 ~ 4 n> 


r.  i.  2 ~*t>  _i_  icu+cji 

L>?  r,  J 


. (40) 


These  two  functions  constitute  the  solution  for  the  case  cf  the  force 


at  (6  0fc  ) normal  to  the  plane  boundary. 


ie  Bounc 


In  this  case  we  take  * Fx  - 0 and  -0  • The  boundary  condi- 
tions then  become,  on  2.-0  , 


^ = _ *g].0  (41) 

= 27T^[V'2,')(§§i  + If5]  - xfitz-  $r]-  ^ (42) 

<s,  - - *l£§  - i&r]  -*  • to) 

Differentiating  (42)  with  respect  toy  and  (43)  with  respect  to  x and  sub- 


tracting, we  find,  on  z-O  , 


or  - 0 


Hence,  on  z -0  , 


9B* 


= £ 


-7- 


Also,  in  2 ? 0 , 


V1  = __  J_  ®J5i 

V “2>z  p 3z 


Honce,  from  (13), 


13l 

3 z 


4tt//  J0-  "1^“  dv 

4^r  ^ ^ x 


by  (32).  Thus 


B 


X 


Prom  (43 ) , on  z = 0 , 


»jL 

32 


(45) 


by  (44).  Also,  in  z ^ O 

v'  - |#-J 


x o Fm 

>/  22 


Hence , by  (13 ) 


(46) 


But  the  right  hand  side  of  (46)  vanishes  since  £ — *0  as  T— *0.  Renee 

(l-2v)Bz  - |f-  =0  (47) 

throughout  the  region  Z 9 0 . 
ft*on  (41 )»  on  2 = 0 , 


-8- 


4-  *)  If 


- x 


(42) 


3 2' 


Bx 


4%  -* 

Z>  Z* 


Now,  on  z = 0 , we  liave,  from  (45), 


2V 


d B, 
3x 


— x 


3XB,  _ 

" 


where  P*  *■  x*  + + C*  . Bat,  on  Z = 0 , 

ft  _ 2v)  ? *0  % x 

('  1V)^5%  27T/4P. 


and 


P,C  9*  / 1 \ = .?f?  cx X 
2ttju  vrtzy I?  ) 27I>/  ^ 


Hence,  (42)  may  be  rewritten  as 


x = 4-v)^-|^_Hv;|S_ 


3R  c.3  x 


.&£_  /j_\.  <, 

27r>£/  axaz^j  u 


(49) 


on  Z - 0 . The  Lap  la  clan  of  the  left  side  of  (49)  is,  in  the  region  Z ^ 0 , 


v1*  - 


X 2>V. 


>/  IF 


I-2V  T>F* 


yU  3K 


Hence , from  (13 ) , 


X = *rhrf&[<S  - 0-2r)-§£-]  dv  . 


(50) 


The  first  term  in  the  integrand  vanishes  since  £ — * 0 c, s T -*  0,  The  second 
term  is  integrated  by  parts  and  the  surface  integral  vanishes,  leaving 


X - i%r  JV, 


3G 

"8T 


d\J 
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which,  by  (17)  and  (25),  is 


O-i v)r  _a_  h L 

4-TT/4  2 X (^F,  f?x 


(51), 


Then,  froo  (51),  (49)  and  (45), 


e^_  a/g  _ __  (i-2v)^x 

2TTJJF?X 


2 Z' 


3Rcx(z  +c) 
1%M 


(52) 


or 


2(i-»h 


df_  = 0-ir)Rx  _ P,  c /, 
2z  ZTyU  rx(kx+*+<)  TW/7W7 


(51) 


finally,  from  (53)  and  (47),  we  hare 

o _ (l-2v)FSx  PmCX 

* nr^Rx(jf+z  +c)  2izm 


(54) 


I* 


(l-2  yT  P*  x , Q-i  v)P,  cx 

1TTJ/(Rx  + z +C)  2.TCM  fe  + z +c) 


(55) 


These  two  functions,  in  addition  to  (45),  coni. rise  the  solution  for  the 
case  of  the  force  at  (CtOt  c ) parallel  to  the  plane  boundary. 


Carper! son  with  Previous  Results 


The  previous  solution,  mentioned  in  the  Introduction,  was  given  in 
terms  of  the  Galerkin  vector  F (not  to  be  confused  with  the  body  force  F 
in  the  present  paper). 

For  the  case  of  a force  normal  to  the  plane  boundary  the  solution 
obtained  was 


-10- 


(56) 


wihrfc  +[*N-']'s 


- 2 c 2/^ 

] 'Uf(*x 


2 ~ VC 


+ z + c 


'} 


and,  for  the  case  of  a force  parallel  to  the  plane  boundary, 

e ‘&fcr  {*.  + f?.-2cA 

+ 4{l-v)(l-Jv)J(z  +c)>yfo  + z +c) 

The  relation  between  the  Galerkin  and  Papkovitch  functions  has  been 
shown  to  be4 

= (MV*F  (58) 

- (|-V)(l7-F  -r-^F)  (59) 

By  Inserting  (56)  and  (57)  in  (58)  and  (59),  it  nay  be  verified  that 
the  previous  and  present  solutions  are  identical. 


-11- 


